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SUMMARY 

First-order  relationships  between  the  position  and  velocity  errors  at 
start  and  end  of  a  Hohmann  transfer  are  derived  in  a  two-dimensional  analysis. 
The  results  are  then  used  to  evaluate  the  effects  of  uncorrected  first  injec¬ 
tion  errors  upon  the  mean  radius  and  eccentricity  of  the  final  (nominally 
circular)  orbit,  assuming  perfect  application  of  the  final  impulsive  velocity. 

Finally  the  mean  radius  and  eccentricity  errors  are  expressed  in  terms 
of  post-orbital  velocities  required  for  nulling,  such  as  would  be  required  in 
station-keeping  communication  satellite  systems. 
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1  INTRODUCTION 

When  placing  satellites  into  high  circular  orbits,  the  flight  plan 
usually  adopted  consists  of  two  phases  of  propelled  motion  (neglecting  any 
specific  requirements  for  parking)#  The  first  phase  extends  from  launch  to  a 
condition  in  which  propulsion  ceases,  resulting  in  the  launch  vehicle  entering 
an  elliptic  orbit  with  a  low  perigee  height  and  apogee  close  to  the  height  of 
the  final  orbit.  In  practice  the  vehicle  will  enter  the  intermediate  or 
"transfer51  orbit  close  to  perigee^  and  will  coast  until,  just  before  apogee 
is  reached,  the  second  propulsive  phase  begins.  This  phase  (of  comparatively 
short  duration  compared  with  that  of  coasting  transfer)  continues  until  the 
vehicle  has  attained  circular  orbital  speed  at  which  time  final  injection  is 
considered  to  have  occurred. 

High  circular  orbits  are  of  prime  interest  in  communication  satellite 
systems,  and  in  such  cases  it  is  often  operationally  difficult  to  monitor  the 
parameters  of  the  transfer  ellipse,  since  this  requires  coverage  by  ground 
tracking  stations*  Consequently  injection  errors  at  the  start  of  transfer 
(velocity  and  position)  may  remain  uncorrected,  resulting  in  final  orbit  errors 
additional  to  those  incurred  by  imperfect  propulsion  during  the  second  phase. 

The  object  of  this  report  is  to  determine  the  effect  of  errors  at  the 

start  of  transfer  upon  the  errors  in  the  major  parameters  of  the  final 

(nominally  circular)  orbit.  The  analysis  is  two-dimensional,  and  consequently 

the  results  are  confined  to  in-plane  errors.  In  order  to  simplify  the  analysis 

considera/tion  is  confined  to  the  special  case  of  Hohmann  transfer.  In  this 
1  2 

special  case  9  coasting  starts  at  perigee  of  the  transfer  ellipse  and  continues 
until  apogee  is  reached,  where  application  of  a  horizontal  velocity  impulse 
completes  the  injection. 

The  analysis  contained  in  this  report  is  divided  into  three  logical 
phases. 

(i)  The  derivation  of  first-order  relationships  between  position  and 

velocity  errors  at  start  of  transfer  to  those  at  the  end  of  transfer,  using 

4,5 

the  well-established  properties  of  ballistic  trajectories  •  In  the  first 
instance  the  relationships  are  expressed  in  a  polar  earth-centred  axis  system, 
but  for  the  sake  of  completeness,  these  results  are  later  transformed  into  a 
cartesian  reference  system.  The  latter  system  is  likely  to  be  the  most  suit¬ 
able  one  in  which  to  evaluate  the  error  contributions  of  an  inertial  guidance 
system. 
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1: 


(ii)  From  the  first  phase  it  is  possible  to  evaluate  the  consequent 
injection  errors  into  the  final  orbit,  assuming  that  the  only  errors  are  those 
at  the  start  of  transfer  (i.e.  the  apogee  impulse  is  applied  as  in  the  case  of 
error-free  transfer).  Final  injection  errors  are  then  transformed  into  the 
equivalent  errors  of  the  mean  radius  and  eccentricity  of  the  final  (nominally 
circular)  orbit. 

(iii)  Since  many  communication  satellite  systems  involve  accurate 

station-keeping  for  periods  of  years,  it  follows  that  such  systems  must  be 

placed  in  orbits  of  exceedingly  high  accuracy.  In  order  to  achieve  this 

accuracy  in  practice  it  vri.ll  be  necessary  to  make  post-orbital  velocity  ad jus t- 
7 

ments  by  means  of  small  thrust  jet  forces,  located  an  the  satellite  itself* 

The  velocity  requirement  from  such  a  system  will  reflect  in  a  direct  loss  of 
potential  payload  for  communication  purposes,  and  it  is  therefore  useful  to 
transform  the  orbit  parameter  errors  of  (ii)  into  equivalent  velocity  correc¬ 
tions.  The  corrections  are  defined  as  those  necessary  to  null  the  mean  radius 
and  eccentricity  errors  of  the  final  orbit. 

2  AS J3F  THE  /Jk\LYS IS 

2  •  1  Axi s^  sypte mis 

As  has  already  been  mentioned,  the  present  analysis  is  limited  to  two 
dimensions,  and  consequently  docs  not  cover  the  evaluation  of  the  effects  of 
errors  normal  to  the  plane  of  injection  (which  is  also  assumed  to  be  the  plane 
containing  the  final  orbit).  The  primary  analysis  is  undertaken  in  a  polar 
axis  system  (see  Fig.1  )  in  which  position  is  defined  by  two  co-ordinates 


and 


r  geocentric  radius 
0  angular  range 


and  velocity  is  defined  by  a  further  two 


and 


V  magnitude  of  velocity  vector 

G  direction  of  vector  relative  to  the  local  horizontal 
(termed  K climb  angle”) 


2.2  The.  I-Iojimann  trjvns f  er 

The  coasting  transfer  is  considered  to  consist  of  an  ellipse  (termed 
“transfer  ellipse”).  The  first  phase  of  propulsion  is  assumed  to  terminate 
(in  the  error-free  or  nominal  case)  at  the  perigee  of  the  ellipse,  and  the 
second  phase  begins  at  the  apogee  of  the  ellipse.  Thus  the  vehicle  coasts 
over  exactly  half  of  the  ellipse.  The  second  propulsive  phase  is  taken  as 
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consisting  of  an  impulsive  addition  of  velocity  (locally  horizontal  in  the 
error-free  case),  such  that  after  the  impulse  has  been  applied  the  vehicle  has 
achieved  circular  orbital  velocity.  At  this  point  the  vehicle  is  considered 

*  to  be  injected  into  the  final  orbit,  and  in  the  error-free  case  this  orbit  will 

1 

be, exactly  circular.  The  transfer  just  described  was  first  proposed  by  Hohmann  . 

1  2 

*  The  mathematical  characteristics  of  such  a  transfer  are  well  known  *  and  are 
listed  in  Appendix  A,  and  illustrated  in  Fig. 2. 

At  first  sight  the  assumption  of  using  a  Hohmann  transfer  in  practice  may 
seem  to  be  rather  sweeping.  However,  many  practical  launching  systems  will 
employ  transfers  which  are  close  approximations  to  this  ideal  (e.g.  the  ELDO 
system  described  in  Ref .3).  Firstly,  for  optimum  performance  with  real  (finite 
thrust)  vehicles  it  transpires  that  entry  into  the  transfer  ellipse  occurs  very 
close  to  perigee,  especially  if  no  altitude  restrictions  are  placed  on  the  first 
injection  point  from  the  point  of  view  of  obtaining  satisfactory  radio  tracking 
coverage.  Secondly,  the  finite  acceleration  encountered  during  the  second 
r  propulsive  phase  is  likely  to  be  sufficiently  high  to  approximate  this  phase  to 

an  impulse,  especially  if  a  spin-stabilised  solid-propellant  stage  is  used. 
Typically  the  angular  range  change  during  propulsion  may  be  5°,  compared  with 
180°  during  the  coasting  transfer. 

In  view  of  these  practical  considerations,  a  Hohmann  transfer  is  taken  as 
the  reference  trajectory  in  the  analysis  of  this  report.  The  resulting  simplifica¬ 
tion  is  considerable  and  the  loss  of  accuracy  in  predicting  the  precise  effect  of 
injection  errors  negligible. 

2.3  Magnitude  of  the  errors 

In  the  whole  analysis  of  this  report  (see  Appendices  A,  B  and  C)  the  rela¬ 
tionships  between  input  disturbances  (i.e.  position  and  velocity  errors  at  start 
of  transfer)  and  output  quantities  (i.e.  errors  at  end  of  transfer  or  final  orbit 
parameter  errors  or  necessary  corrective  velocities)  are  linearised.  Therefore 
the  resulting  equations  can  only  be  applied  strictly  if  the  perturbations  are 
v  infinitely  small,  since  second  and  higher  order  terms  are  neglected  by  implica¬ 

tion.  The  approximation  involved  by  such  procedures  is  not  likely  to  be  serious 
»  however,  since  in  practical  launching  systems  the  velocity  injection  errors  at  the 

start  of  transfer  are  likely  to  be  of  the  order  of  tens  of  ft/sec  compared  with 
actual  velocities  of  several  thousands. 

2.4  Application,  of  impulsive  velocity  during  the  second  burning  phase 

In  the  analysis  it  is  assumed  that  the  velocity  impulse  is  added  perfectly 
i.e.  we  are  only  investigating  the  effect  of  first  injection  errors  upon  the 
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final  orbit*  However,  it  is  necessary  to  clarify  what  is  meant  by  "perfect1* 

application  -  of  the  impulse*  In  the  nominal  case  (i.e.  no  injection  errors  at 

the  start  of  transfer)  the  impulse  is  applied  at  a  time  Tq/2  after  start  of 

transfer  (where  T^  is  the  orbital  period  of  the  nominal  transfer  ellipse)  in 

a  direction  which  is  simultaneously  parallel  to  the  local  horizontal  and  normal 

to  the  radius  vector  at  the  start  of  transfer  (see  Fig.2).  The  magnitude  is 

(V  -V  )  where  V  is  the  circular  orbital  speed  at  the  final  radius  r  and  V 
°2  ^  °2  2  2 

is  the  velocity  at  the  end  (apogee)  of  transfer  (all  values  nominal).  In  the 
general  perturbed  case  however,  i.e.  when  first  ejection  errors  are  present, 
the  local  horizontal  at  the  end  of  transfer  (transfer  being  still  considered  to 
terminate  after  a  time  T^/2  since  no  more  tracking  information  is  available) 
does  not  coincide  with  the  direction  nominal  to  the  radius  vector  at  nominal 
transfer  start  (e.g.  see  Fig. A),  due  to  position  errors  at  the  end  of  transfer* 
We  must  therefore  differentiate  between  two  cases  of  "perfect**  impulse  applica¬ 
tion:-  the  first  one  in  which  the  impulse  is  aligned  parallel  to  the  local 
horizontal  at  time  T  /2  after  start  of  transfer,  and  the  second  one  in  which 
it  is  aligned  in  a  space-fixed  direction  normal  to  the  radius  vector  at  nominal 
transfer  start.  Both  methods  ox  alignment  are  covered  in  this  analysis  since 
both  can  be  achieved  in  practice.  For  example  the  first  one  is  appropriate  to 
alignment  using  sensors  detecting  the  local  horizontal  (optical  horizon  sensors), 
and  the  second  method  would  be  realised  using  either  inertial  guidance  or  a  spin 
stabilised  final  burning  stage. 

3  RESULTS, 

3.1  position  jand_  velocity  errors  at  start _of  transfer 

to  jjillPJL0  at  the  end 

3*1.1  Polar  axi^sy^em 

The  detailed  analysis  is  contained  in  Appendix  A.  The  results  enable  the 
perturbations  at  beginning  and  end  of  transfer  to  be  related  in  the  form 

I  &y 2 1  =  MlSyj  0 ) 

where  ] 6y |  is  the  perturbation  state  vector  at  the  start  and  | 5y^ |  that  at  the 
end.  Each  of  these  vectors  is  f our-dimensional,  containing  two  position 
co-ordinates  (r,  $>)  and  two  velocity  co-ordinates  (V,  6)  as  shown  in  Fig. 2, 

The  4  x  4  matrix  |m|  contains  sixteen  derivatives,  all  of  which  are  evaluated 
in  Appendix  A  and  are  tabulated,  for  the  sake  of  completeness,  in  Table  l(a). 

The  derivatives  (other  than  those  which  are  scro  or  unity)  are  plotted  in 
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Figs, 11,  12  and  13.  The  evaluation  is  in  terms  of  a  dimensionless  parameter, 
termed  "transfer  ratio11 


(2) 


where  r .  is  the  radius  at  start  of  transfer 
1 

r^  is  the  radius  at  end  of  transfer  (also  equal  to  the  radius  of  the 
final  circular  orbit)# 

For  generality  to  derivatives  and  normalised  with  respect  to  quantities  which 

are  dependent  only  upon  the  initial  radius  e#g#  r^  and  where 

1 

V2  =  Ji-  (3) 

0'] 

Vq  is  the  circular  orbital  speed  at  r. 

and  p  is  the  earth’s  gravitational  constant. 

In  Figs.llj  12  and  13  the  derivatives  are  plotted  non-dimensionally,  but 
additional  dimensional  scales  are  provided  for  the  special  case  in  which,  the 
altitude  at  start  of  transfer  is  300  nm.  The  numerical  values  taken  in  deriving 
the  scaled  derivatives  are 


Earth  radius  3437*75  nautical  miles 
ji  =  1.40673  x  1016  ff/sec2 
1  nm  =  6080  ft  . 

3.1.2  Cartesian  axis  system 

The  results  obtained  in  Appendix  A  have  been  transformed  into  an  cartesian 
axis  system,  as  shown  in  Fig. 2.  The  detailed  analysis  of  the  transformation  is 
contained  in  Appendix  B  and  the  results  summarised  in  Table  l(b).  The  geometrical 
relationships  between  small  perturbations  in  polar  and  cartesian  axis  are 
illustrated  in  Fig. 8. 

3*2  Derivatives  relating  position  and  velocity  errors  at  start  of  transfer 

to  mean  radius  and  eccentricity  errors  of  final  orbit 

In  Appendix  C,  errors  at  the  end  of  transfer  are  related  to  injection 
errors  for  the  final  orbit  (assuming  the  two  cases  of  "perfect"  impulse 
previously  described)  and  hence  to  errors  in  mean  radius  and  eccentricity. 
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The  results  are  shown  (both  dimensionally  and  non-dimensionally)  in 
Figs. 14,  11)  11  and  20.  Since,  as  is  shown  in  Appendix  C,  eccentricity  is 

a  non-linear  function  of  injection  errors,  we  cannot  apply  the  usual  formula 


8e 


( j3e 

a  $6 


Ha 


(4) 


In  the  presence  of  a  single  input  error  it  can  be  shown  that  a  linear  relation¬ 
ship  will  exist  between  input  err* or  and  resulting  eccentricity,  provided  the 
input  error  is  replaced  by  its  modulus,  since  positive  eccentricity  is  produced 
by  either  positive  or  negative  inputs,  i#e# 


6e  n 


de  \ 


dV 


1/ 


6V  etc. 
1 1 


(5) 


where  (de/dV,|)  is  positive. 

In  the  presence  of  multiple  input  errors,  equation  (4)  must  be  replaced 
by  a  more  complex  expression  which  will  not  be  pursued  here# 

3  #3  Derivatives  relating  position  and  velocity  errors  at  s  tarjp  o_f_  transfer  • 

velocity  corrections  required  to  null  orbital  errors 

In  Appendix  C  the  relations  between  post-orbital  corrective  velocity 

requirements  and  orbital  errors  are  established,  thereby  making  it  possible  to 

derive  expressions  for  the  corrective  velocities  in  terms  of  errors  at  the 

start  of  transfer#  Two  velocity  corrections  arc  calculated,  6u  that  necessary 

a 

to  null  mean  radius  error,  and  5ue>  that  for  correcting  eccentricity#  Further¬ 
more  it  is  shown  that  both  radius  and  eccentricity  errors  can  be  eliminated  by 
a  total  impulse  application  equal  to  either  Su^  or  8u^,  which  ever  is  the 
larger  numerical  quantity* 

Once  again,  since  5u  is  a  non-linear  function  of  6  V, ,  80  ,  Sr.  and  80 

e  i  i  l  i 

(even  fox'  small  perturbations)  it  cannot  be  evaluated  from  a  linear  relation 

of  the  form  of  equation  (a)  when  multiple  inputs  are  present  at  the  start  of  , 

transfer . 

Corrective  velocity  derivatives  are  plotted  in  Figs. 1 6,  18  and  23.  In  « 

Fig.18,  du  /dd  is  normalised  by  dividing  by  V.  (rather  than  V  )  since  V. 80. 

e  i  1  o^  1  1 

is  directly  vertical  velocity  error.  Fence  Figs. 1 6  and  18  now  form  a  direct 

comparison  of  the  relative  effects  of  vertical  and  horizontal  cut-off  velocity 

errors.  Similarly,  in  Fig. 23  the  scaled  quantity  plotted  in  the  ordinate  is 

(due/r_j  .  Nov/  r^6^  is  the  longitudinal  position  error,  so  that  again  we 


can  compare  the  relative  effects  of  position  errors  in  two  directions  from 
Figs. 21  and  23. 

4-  DISCUSSION 

4 * 1  Relations  between  errors  at  beginning  and  end  of  transfer 

There  is  little  comment  that  can  be  passed  upon  these  results,  except  to 
mention  their  possible  use  in  the  overall  detailed  evaluation  of  any  launcher 
guidance  system  that  contains  a  Hohmann  transfer  (or  a  close  approximation)  in 
the  flight  plan.  It  was  for  this  type  of  evaluation  that  the  derivatives  in 
Appendix  A  were  converted  into  a  cartesian  axis  system.  Such  a  system  of  axes 
is  the  obvious  choice  for  detailed  investigation  of  inertial  guidance  errors, 
since  most  inertial  navigators  for  space  launchers  will  measure  and  navigate  in 
these  axes  (e.g.  Section  8  of  Ref. 7). 

It  is  worth  noting  that  initial  errors  in  50^  and  8$^  perturb  the  transfer 
ellipse  in  a  different  way  to  errors  in  8V^  and  8r^ •  The  former  produce  no 
first  order  change  in  the  mean  radius  of  that  orbit  (resulting  in  zero  change 
of  radius  or  velocity  at  the  end  of  transfer)  while  the  latter  produce  changes 
in  most  of  the  transfer  parameters  (resulting  in  changes  to  all  the  co-ordinates 
at  the  end  of  transfer). 

4.2  Errors  in  parameters  of  the  final  orbit 

Errors  in  mean  radius  (6a)  are  only  produced  by  5Vj  and  5r^  inputs 
(Figs, 14  and  19),  and  these  errors  are  independent  of  the  method  of  aligning 
the  apogee  impulse  (this  only  affects  vertical  velocity  errors  at  injection). 

The  mean  height  error  for  a  given  initial  velocity  error  increases  with  transfer 
ratio,  as  does  mean  height  error  for  a  given  initial  radius  error.  Typical 
values  for  injection  into  a  geo-stationary  orbit  (19450  nm  altitude),  assuming 
transfer  to  start  at  300  nm  altitude  are 


Eccentricity  8e  results  from  any  type  of  initial  perturbation  (Figs. 15, 

17,  20  and  22).  For  5Vj  and  8r^ ,  5e,  for  a  given  initial  disturbance,  increases 
with  transfer  ratio.  Also  eccentricity  can  be  minimised  by  aligning  the  apogee 
impulse  to  a  fixed  space  direction  rather  than  to  the  local  horizontal,  since 
in  the  former  case  the  vertical  velocity  error  at  the  end  of  transfer  is  partially 
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compensated  by  a  locally  opposite  application  of  velocity  during  the  impulse 
(i.e.  is  of  opposite  sign  to  and  dd^/dr^  is  of  opposite  sign 

to  The  effect  of  60 ^  or  8<j^  errors  is  considerably  different  from 

that  of  6V  or  8r^ .  Firstly,  the  increasing  effect  with  increasing  transfer 
ratio  is  much  less  marked  ( [de/80^ ]^_  actually  decreases  with  increasing  n, *) 
and  secondly,  the  minimum  eccentricity  error  is  produced  by  aligning  the 
apogee  impulse  to  the  local  horizontal,  rather  than  to  a  fixed  space  direction 
(actually  (de/d^)  q  is  zero  for  all  value  of  n). 

For  the  geo-stationary  orbit  eccentricity  derivatives  are 


=  O.OOO56  per  ft/sec 

Vk=o 


0e 


av. 


1/k=1 


=  0.00037  per  ft/sec 


=  0.00285  per  nm 

r1/k=o 


de\ 

sVk=i 


=  0.0018  per  nm 


ae 


0.09  per  radian  (or  0.0000027  per  ft/sec  of  vertical  velocity  error) 


1/ lc=o 


de 

dfT 


1/k=1 


=  1.17  per  radian  (or  O.OOOO36  per  ft/sec  of  vertical  velocity  error) 


/de\ 

\  ^1/k=1 


0  (or  zero  per  nm  of  range  error) 


0.46  per  radian  (or  0,00012  per  nm  of  range  error) 


We  can  see  from  these  figures  that  horizontal  velocity  error  (6V^)  has  a 
larger  effect  upon  final  orbit  eccentricity  than  vertical  velocity  error 
(Vj  60  )  even  for  k=1  (apogee  impui.se  aligned  in  space),  which  is  the  case  most 
favourable  to  5V^  and  least  favourable  to  60  . 


*The  suffix  lc=o  refers  to  the  system  with  the  apogee  impulse  aligned  to 
the  local  horizontal,  while  k=1  corresponds  to  space  alignment  (see  Appendix  C, 
para.  A3.2)* 
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Similarly  radial  position  error  (8r^ )  has  a  larger  effect  than  range 
position  error  (r^  8^  ) . 

4.3  Post-orbital  corrective  velocities 

Before  discussing  the  numerical  values  (Figs.l6,  18,  21  and  23)  it  is 
appropriate  to  review  the  possible  implications  of  corrective  velocities  upon 
a  practical  communication  satellite  launching  system. 

Firstly  it  must  be  noted  that  orbital  period  must  be  corrected  to  a  very 
high  degree  of  accuracy,  since  relative  station  keeping  between  individual 
satellites  in  the  complete  system  must  be  maintained  to  within  a  few  degrees 
over  a  period  of  several  years.  Since  orbital  period  is  directly  related  to 
mean  radius  (e.g.  Ref .5),  this  implies  that  any  initial  mean  radius  error  (8a) 
must  be  subsequently  corrected  to  an  extremely  high  degree  of  accuracy.  In 
fact  the  final  accuracy  must  be  of  the  order  of  a  fraction  of  an  inch  per  second 
(in  terms  of  velocity  error)  so  that  8u&  is,  in  practice,  a  good  measure  of  the 
contribution  that  can  be  expected  from  first  injection  errors  (other  contribu¬ 
tions  can  be  expected  from  errors  in  the  impulse  itself,  but  these  are  not  the 

subject  of  the  present  investigation).  It  must  be  remembered  that  8u  is 

a 

calculated  here  on  the  basis  of  tangential  impulse,  so  that  it  is  essential 
that  practical  post-orbital  corrections  should  be  capable  of  representation  by 
such  impulses. 

Regarding  correction  of  initial  eccentricity  error  (8e)  the  situation  is 
by  no  means  so  well  defined.  Residual  eccentricity  errors  result  only  in  cyclic 
station  keeping  errors  (an  eccentricity  5e  produces  a  sinusoidal  phase  error  of 
amplitude  i2Se  and  period  equal  to  that  of  the  nominal  orbit).  Some  degree  of 
cyclic  error  may  be  permissible  in  an  actual  communication  system,  although  this 
will  depend  upon  such  details  as  the  permissible  movement  of  the  ground  aerials. 
Therefore  the  values  of  Su^  given  here  represent  an  upper  limit  to  the  corrective 
velocity  that  would  be  required. 

Turning  now  to  the  present  results,  let  us  examine  8u  and  8u  required 
from  an  initial  error  &V^  (Fig.16).  Generally  speaking  both  du^/dV.  and  du 
increase  with  increasing  transfer  ratio,  but  with  some  slight  variation  to  this 
rule  for  n  between  1  and  2.  For  a  synchronous  orbit  8u  is  larger  than  8u  when 

Q,  C 

the  apogee  impulse  is  aligned  in  space,  while  the  reverse  is  true  if  the  impulse 

is  aligned  to  the  local  horizontal#  As  is  explained  in  Appendix  C,  6u  >  6u 

a  0 

implies  that  the  perturbed  final  orbit  does  not  intercept  the  desired  nominal 

circular  orbit,  while  6u  <  6u  implies  that  interception  does  occur#  Further- 

a  e 

more,  if  8u&  >  8ue  it  means  that  eccentricity  error  can  be  nulled  during  the 
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mandatary  process  of  orbital  period  correction.  Also  the  positive  margin 

between  8u  and  6u  allows  some  additional  vertical  velocity  error  to  occur 
a  e 

(e.g.  from  mis -direction  of  the  apogee  impulse)  and  still  permit  eccentricity 
nulling  within  the  capability  of  Su  . 

~  cl 

From  Fig.l6,  du^/<3V^  =  2.2  in  a  synchronous  orbit*  This  emphasizes  the 
importance  of  keeping  uncorrected  first  injection  forward  velocity  errors  as 
small  as  possible,  since  their  effect  is  magnified  in  relation  to  post-injection 
velocity  requirements* 

For  an  initial  vertical  velocity  error  (V^  60^ )  (see  Fig.18),  no  correction 

is  necessary  for  mean  radius  (bu^  =  0),  and  the  correction  for  eccentricity  is 

quite  small  (l/V.  du  /60  =  0.18  in  the  most  unfavourable  impulse  alignment 

i  e  I 

case)* 

For  an  initial  radius  error  (6r^)  (see  Fig#2l),  the  corrective  velocities 
are  similar  in  form  to  those  for  6Vj  •  In  a  synchronous  orbit 

du 

=  11.5  ft/sec  per  nm 


and  the  eccentricity  correction  is  slightly  larger  for  the  case  of  horizontal 
impulse  alignment,  and  smaller  for  space  alignment. 

For  a  range  error  (r^ §9^  )  there  is  no  period  correction  required,  and  no 
eccentricity  correction  in  the  case  of  horizontal  impulse  alignment.  For  the 
case  of  space  alignment  du^/r^d^  =  0*62  ft/sec  per  nm  (synchronous  orbit) 
indicating  that  the  effect  of  range  errors  is  considerably  smaller  than  that 
of  height  error. 


Orbital  errors  and  corrective  velocities  in  the  presence  of 


In  any  practical  launching  system  errors  will  exist  in  all  four  co-ordinates 
at  the  first  injection  point  (in  addition  to  errors  of  the  apogee  impulse). 

Usually  these  individual  errors  will  be  defined  statistically  and  it  is  desired 
to  evaluate  the  effect  of  such  errors  upon  final  orbit  parameters  and  corrective 
velocities.  Such  effects  will  again  be  derived  in  statistical  form*  Often  it 
is  possible  to  suppose  that  the  "input"  quantities  are  normally  distributed,  and 
uncorrelatcd  with  one  another*  If  the  "input"/" output"  relationships  are  linear 
it  is  easy  to  show  that  the  "output"  quantities  are  also  normally  distributed  and 
the  standard  deviation  of  the  "output"  is  readily  obtained  in  terms  of  the  "input" 
standard  deviation  and  the  appropriate  linear  relationship.  If  the  relationship 
is  non-linear  the  statistical  evaluation  of  "output"  is  more  complex  and  no 
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general  solution  exists.  Unfortunately,  such  is  the  case  in  the  case  of  our 
"outputs",  eccentricity  (5e)  and  corrective  velocity  (5u  )•  This  problem  has 
already  been  encountered  in  Ref. 7.  No  solution  is  suggested  here,  the  sole 
object  of  mentioning  it  being  to  advise  exercise  of  caution  in  applying  the 
results  of  this  report. 

5  CONCLUSIONS 

(1)  Derivatives  have  been  evaluated,  relating  errors  (position  and 
velocity)  at  the  start  of  a  Hohraann  transfer  to  those  at  the  end.  Such 
derivatives  are  of  use  in  the  detailed  investigation  of  the  accuracy  of  any 
guidance  system  employed  in  injecting  satellites  into  high  circular  orbits 
(such  as  those  required  for  communication  satellite  systems). 

(2)  These  derivatives  have  been  used  to  determine  the  effect  of  single 
errors  at  first  injection  upon  the  parameters  of  the  final  orbit,  and  the 
resulting  necessary  post  injection  corrective  velocities.  For  large  values  of 
transfer  ratio  (final  orbit  radius  divided  by  first  injection  radius)  it  has 
been  shown  that  horizontal  velocity  errors  are  considerably  more  significant 
than  vertical  ones,  while  radial  position  errors  have  a  larger  effect  than 
range  position  errors.  As  an  example,  for  injection  into  a  geo-stationary 
orbit  each  unit  of  horizontal  velocity  error  will  require  at  least  2.2  units 
of  post  orbital  correction  and  each  nautical  mile  of  radial  position  error 
needs  11.5  ft/sec  correction  after  injection. 
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Appendix  A 

PROPERTIES  OF  THE  HOHMANN  TRANSFER,  INCLUDING-  RELATIONS  BETWEEN 
SMALL  PERTURBATIONS  AT  THE  BE&DRJIHO  AM)  END  OF  THE 
TRANSFER  PHASE,  DERIVED  IN  A  POLAR  AXIS  SYSTEM 

A • 1  General  properties  of_  the  Eohroann  transfer  orbit 

The  ideal  transfer  (described  in  Ref 3.1 ,  2)  consists  of  ballistic  flight, 
extending  over  half  of  an  ellipse,  commencing  at  perigee  and  ending  at  apogee 
(Fig. 2),  The  initial  and  final  velocities  of  the  transfer  are  defined  in  terms 
of  the  earth’s  gravitational  constant  (p)  and  the  initial  and  final  geocentric 
radii  (r^  and  ) .  TJie  appropriate  relations  are 


V,  =  V  r,  =  V 


2n 


0^  n+1 


(A.1) 


Y2  -  \P2 


V 


!  2 

°2  i  'n+T 


(A#2 ) 


where  Vq  =  Vp/r^  is  the  circular  orbital  speed  at  the  initial  radius 

and  V  =  Vu/ r_  is  the  circular  orbital  speed  at  the  final  radius  zr  • 

o2  "  2  ^  2 

The  parameter  p  has  been  previously  used  in  B.ef#4  and  is  defined  as  the  square 
of  the  ratio  of  speed  to  the  circular  orbital  speed  at  the  appropriate  radius* 


It  will  be  noted  that 


Ii 

v„ 


v2 


V 


=  n 


where  n  =  r„/r, . 

2  1 

The  dimensionless  parameter  n  (ratio  of  final  to  initial  radii)  will  be 
termed  the  transfer  ratio,  and  all  derivatives  relating  initial  and  final 
perturbations  will  be  expressed  as  a  function  of  this  quantity. 

The  transfer  time  is  given  by  Tq/2,  Tq  being  the  orbital  period  of  the 

transfer  ellipse.  Now  orbital  period  is  related  simply  to  the  mean  radius  of 

5 

the  orbit 


T  =  2x  1— 
o  -J  ( 1 


(A.3) 


a  =  i(r1 


r2)  • 


(A.4) 


where 
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The  establishment  of  a  circular  orbit  at  radius  r^  is  effected  by  adding 

a  locally  horizontal  velocity  impulse  (of  magnitude  V  -V  )  on  arrival  at  the 

°2  2 

end  of  the  transfer# 

A»2  The  effect  of  a  small  velocity  perturbation  6V^  at  the  start  of  transfer 

Consider  the  situation  presented  in  Fig. 4#  Point  S  represents  the  start 

of  the  transfer  and  point  E  is  the  end  point  of  the  nominal  (i#e.  unperturbed) 

transfer.  Now  point  A  represents  the  apogee  of  the  perturbed  transfer  ellipse# 

As  illustrated,  A  is  shown  for  a  positive  perturbation  in  •  Since  A  will 

occur  later  in  time  than  E,  the  actual  perturbed  transfer  will  end  at  point  E', 

assuming  the  impulsive  velocity  (V  -V  )  is  added  at  the  instant  corresponding 

°2  2 

to  the  end  of  the  nominal  transfer. 

A. 2.1  V  and  r  perturbation  at  end  of  transfer 

Let  5t  be  the  perturbation  in  the  half -period  of  the  transfer  ellipse,  due 

to  initial  perturbation  6V. •  Since  r  =  0  at  point  A,  the  difference  of  radius 

*  *.2 

vector  between  A  and  E'  is  of  order  r(St)  ,  where  r  is  described  near  to  apogee. 
Since  6t  will  be  proportional  to  8V,j  (for  small  6V^),  we  can  neglect  the  differ¬ 
ence  between  radius  vectors  at  A  and  Er  when  evaluating  derivatives  describing 
the  first  order  relationship  between  8V^  and  radius  perturbation  at  E*.  Similarly, 
since  radii  changes  between  A  and  E*  are  of  2nd  order,  compared  with  radii  changes 
between  E  and  A,  we  may  also  conclude  that  to  a  first  order  the  velocity  magnitude 
perturbation  at  E^SVg)  equals  that  at  A.  This  follows  from  the  conservation  of 
total  energy  (kinetic  plus  potential)  between  A  and  E*. 

T7e  may  therefore  conclude  that  the  velocity  and  radius  perturbations  at 

E*  are  equal  to  those  at  A,  provided  that  8V,  is  small.  Let  us  assume  that  the 

perturbations  at  A  or  E*  are  5V  ,  5r_  respectively.  Then  from  the  conservation 

22  6 

of  angular  momentum  between  perigee  and  apogee  ,  we  have 
V^r.  =  ^2r2  (norQinal  case) 

(V,j+6Vj)r.|  =  (V2+  8V2)  (r2+ Sr2)  (perturbed  case)  . 

By  differencing  these  equations:- 

r16V1  =  r26V2  +  .  (A.5) 

Note  that  again  we  have  neglected  the  second  order  term  6V^ 
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1 


Similarly  from  the  conservation  of  total  energy  between  apogee  and 
perigee. 


-2-  V?  -  £• 


-2-(V,  +  SV  r  -  J±- 
1  r^ 


2  V2 


V  -  (nominal  case) 


2  n 

=  ^'(Y2+S'VV  "  "(7" 5r”7  (Perturbed  case) 


Again,  by  differencing  and  neglecting  2nd  order  terms 


VV  =  V  8V  +  6r2 

r2 

By  solving  equations  (5)  and  (6)  simultaneously. 


6V_  =  5V, 


and 


5r, 


6V. 


/V  r  V  r  r  • 

/T2  _  _ 2 

V  H  H 


1  - 


V2  r 
2  2 


(A.  6) 


(A. 7) 


(A.8) 


The  derivatives  resulting  from  these  two  equations  can  be  written  in  a 
more  convenient ,  and  non-dimensional  form,  using  the  formula  for  start  and  end 
velocities  as  given  in  equations  (l )  and  (2):- 


av2 

ew 


- .  1 


=  -  2  + 


n 


(A.9) 


V 


dr2 

ev. 


2n(l+n)" 


(A.IO) 


A. 2*2  6  and  perturbation  at  the  end  of  transfer 

Y!e  have  already  seen  that  the  apogee  of  the  perturbed  transfer  ellipse  is 
attained  at  a  later  time  than  the  end  of  the  nominal  transfer  (for  positive  ) , 
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This  is  because  the  mean  height  of  the  perturbed  transfer  ellipse  is  higher 

than  that  of  the  nominal  one.  If  the  perturbed  apogee  occurs  5t  later  than 

i  5T 

that  for  the  nominal,  then  5t  =  2  0,  where  8Tq  is  the  increase  in  orbital 

period  of  the  perturbed  ellipse.  Using  equations  (3)  and  (4)  it  is  possible 
to  relate  5t  to  8r^  and  hence  to  8Vj  (for  small  perturbations).  The  result 
is 


5t 


(A.11) 


Now  one  of  the  equations  of  motion  describing  ballistic  flight  is 


0 

ve 


v2 


*  U, 

=  ~  o  8  + 

r2 


COS  0 


Also  0  =  0  at  point  A,  so  that  at  E* 


(A.12) 


Combining  equations  (ll)  and  (12)  we  obtain  the  following  non-dimensional 
derivative  relating  terminal  climb  angle  error  to  initial  velocity  error 


(A.13) 


To  obtain  6$^  at  the  end  of  the  perturbed  transfer,  we  multiply  by  8t  (the 
negative  sign  indicating  that  6$  is  negative  for  positive  8V. ).  Again  we 

f  0  ,  s  ' 

obtain  a  normalised  derivative  (using  =  V^/r^  ‘ *” 


(A.14) 


The  derivatives  given  by  equations  (9),  (lO),  (13)  and  (14)  are  plotted  in 
Fig.11. 
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A , 3  The  effect  of  a  small  climb  path  perturbation  at  the  start  of  transfer 

This  situation  is  presented  in  Fig.5 •  Point  S  represents  the  start  of 
nominal  and  perturbed  transfers.  E  is  the  end  of  the  nominal  transfer  and  P 
is  the  perigee  of  the  perturbed  transfer.  For  a  positive  60^  perturbation,  P 
will  occur  before  the  start  of  transfer  (as  illustrated  in  Fig.5) •  Let  P  occur 
at  a  time  -8t  (i.e.  6t  before  start,  where  6t  is  not  related  to  the  St  used  in 
para.  A, 2).  Nov;  from  the  equation  of  motion  used  previously,  i.e. 


°  a  V 

VO  =  -  cos  0  +  —  cos  0 

2  r 

r 


o  2  2 

0  near  perigee  equals  (V^/r  -  p/r  )  l/V  ,  which  is  positive,  so  that 


60, 


6t  = 


near  perigee 


X  a  6  < 


6t  = 


50.  V„ 
1  1 


(A.15) 


1  .  JL 
?,  2 
1  rr 

Now  the  radius  r  at  perigee  P  is  of  the  form 


2 

i*e.  of  the  form  r^  -  K!S0  ,  so  that  to  a  first  order  the  perigee  radius  is 
unaltered  from  the  nominal  value  r^ •  Also  the  total  energy  of  the  perturbed 
transfer  orbit  is  unaltered  from  that  of  the  nominal  one  since  the  potential 
energy  is  unaltered  and  the  velocity  vector  at  S  is  only  changed  through  a  small 
angle  80^ •  It  also  follows  that  the  orbital  period  is  unaltered  as  is  the 
apogee  radius  (to  a  first  order)* 

Let  A  be  the  apogee  of  the  perturbed  transfer  (at  radius  r^)*  ^kis 
occurs  at  time  (T  /2  -  ot)  after  start  of  transfer*  Point  E*  represents  the 
end  of  perturbed  transfer,  where  ET  occurs  6t  after  apogee  A*  (As  previously, 
we  are  assuming  that  the  perturbed  transfer  terminates  at  the  same  time  instant 

as  for  the  nominal  transfer,  i.e*  1  /2  after  start*) 

9  o 

A •  3 • 1  V  and  r  perturbation  at  end  of  transfer 

Since  the  total  energy,  perigee  and  apogc-e  radii  are  identical  for  nominal 
and  perturbed  transfers,  V  and  r  at  point  A  will  be  identical  to  the  values  at  E# 
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Also,  since  E*  occurs  at  a  time  6t  after  A,  the  radius  at  E*  will  differ  from 

2  2 

that  A  by  a  term  proportional  to  6t  ,  i.e.  proportional  to  86^ •  Hence  there 
is  no  r  perturbation  at  E’  compared  with  E,  to  a  first  order.  Also,  since 
there  is  no  change  in  total  energy  due  to  perturbation  there  is  no  V  perturba¬ 
tion  at  E*  compared  with  E.  We  can  simply  write 


V 


5 

ee. 


=  o 


(A.16) 


o 


(A.17) 


A»3«2  6  and  $  perturbations  at  end  of  transfer 

Near  apogee,  from  the  previously  used  equation  of  motion 


o 

6 


so  that  at  E’,  80 ^  equals  G  St. 


Erom  equation  (A. 15)  this  gives 


Using  the  previously  established  relations  between  ,  V^,  r^  and  r^  this 
reduces  to  the  dimensionless  form:- 


r 


(A. 18) 


Referring  to  Fig. 5,  the  angular  range  perturbation  at  E*  is  8 0_,  where 
8^  =  -8a  +  8(3. 
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Now 

8t  V1 

6a  =  6t  x  angular  velocity  near  perigee  =  — - — 

r1 

6t  V? 

and  5(3  ~  6t  x  angular  velocity  near  apogee  =  — — — 

r2 

So  that 


(A. 19) 


Substituting  for  6t  from  (A.15) >  we  can  obtain  the  dimensionless  form 


d<t>r 

L 

w. 


2  1  + 


n 


(A.20) 


The  derivatives  and  3^/60^  are  plotted  in  Fig. 12. 

A.4  The  effect  of  a  small  radial  perturbation  at  the  start  of  transfer 

Pig.6  illustrates  the  geometry  of  the  nominal  and  perturbed  transfer;  S 
is  the  nominal  start  point  and  S*  the  perturbed  one.  E  represents  the  nominal 
transfer  end  and  A  the  apogee  of  perturbed  transfer  ellipse.  Since  A  occurs  at 
time  (T^/2  +  6t),  once  again  we  have  the  perturbed  end  occurring  at  El,  some 
distance  before  A# 

As  was  shown  in  para .  A.  2,  the  V  and  r  perturbations  that  exist  at  A  are, 
to  a  first  order,  identical  to  those  at  Ef,  since  Ef  is  close  to  apogee. 

A .4*1  V  and  r  perturbations  at  end  of^tr ansf e r 

Let  the  perturbations  at  A  (or  Ef)  be  b\7^  in  velocity  and  8r^  in  radius. 
Pursuing  the  same  technique  as  was  used  in  para.  A. 2.1  we  have,  from  the 
conservation  of  angular  momentum,  and  total  energy 


V,  r,  =  V  r 
11  2  2 


V  (r  +  8r  )  =  (V  +  8V_)  (r  +  5r  ) 

11  V  2  2'  v  2  2/ 


V,  8r 
1  1 


giving 


r  6V  +  V  8r 
2  2  2  2 


(A. 21 ) 
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and  also 

i  vf  -  ^  =  i  vf  -  n/r2 

£  V2  -  p/r1  +6^  =  1-  (V2  +  6V2)2  -  p/(r2  +  8r2) 

giving 

P-A2  8r1  =  V2  8Y2  +  nAg  5r2  •  (A.  22) 

Solving  equations  (A,2l)  and  (A.22)  simultaneously  we  obtain  the  dimensionless 
derivatives 


so  that 
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This  can  be  reduced  to  the  non-dimensional  form: 


r  . 

!  ri  Tr’  = 


8n 


Vn 


(A. 26) 


Also,  the  perturbed  transfer  end  occurs  at  an  angular  range  with 
perturbation 


-St  x  angular  velocity  near  apogee 


But  the  angular  velocity  near  apogee 


is  V0/r^  so  that 


(A.27) 


Substituting  (A*25)  in  (A«27)>  we  obtain  the  non-dimensional  derivative 


| 

!  r1 


3? 

K 


(_1+n)2  j 
n  Vn  j 


(A. 28) 


The  derivatives  from  equations  (A.23),  (A.24),  (A.26)  and  (A»28)  are  plotted  in 
Fig.13. 


A. 5  The  effect  of  a  small  angular  range  perturbation  at  the  start  of  transfer 

The  effects  of  this  initial  perturbation  are  simple  to  visualise  (see 
Fig.7).  Since  the  apogee  and  perigee  radii  and  the  corresponding  velocities 
are  independent  of  the  angular  position  at  first  injection  it  is  reasonable  to 
conclude  that  the  only  perturbation  to  the  transfer  ellipse  is  in  respect  of 
the  orientation  of  the  major  axis*  This  resulting  change  in  orientation  equals 
the  initial  perturbation  09^  • 


Thus,  quite  simply 


!  j... 

0V 

2 

=  0 

!  V 

1  0 

1 

'd<Aj 

|  A. 

dr 

2 

=  0 

.  ri 

d(f>1 

(a. 29) 


(A.30) 
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ZJ>2 


(A.31) 


(A.32) 
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DERIVATIVES  RELATING-  PERTURBATIONS  AT  START  OF  TRANSFER  TO  THOSE 
AT  THE  END:  CONVERSION  EROI.I  POLAR  TO  CARTESIAN  AXES  SYSTEMS 


B •  1  Definition  of  axes  systems 

The  relevant  systems  are  illustrated  in  Fig.1 .  In  the  polar  system  posi¬ 
tion  is  defined  r  (geocentric  radius)  and  angular  range  $  (the  datum  of  which 
is  arbitrarily  taken  as  the  vertical  at  the  point  of  commencing  the  nominal 
Hohmann  transfer)*  Velocity  is  defined  by  magnitude  V  and  climb  path  relative 
to  local  horizontal,  6* 

In  the  cartesian  system,  position  is  defined  by  two  co-ordinates  x,  and  z* 
These  axes  are  orthogonal  and  inertially  fixed*  The  origin  is  at  the  earth 
centre  0  and  Ox  is  parallel  to  the  velocity  vector  at  nominal  transfer  start* 
Consequently  Os  is  parallel  to  the  <f>  datum  in  the  polar  system  and  is  arbitrarily 
considered  to  be  positive  when  vertically  upwards  at  this  point.  The  nominal 
start  and  end  conditions  for  transfer,  in  the  two  systems  are  obvious,  and  are 
illustrated  in  Fig. 2  for  the  polar  case*  In  cartesian  axes  the  conditions  are 


in  Fig*3>  and  are  as  follows :- 

Initial  x  = 

X1  = 

0 

2  — 

z.  = 

1 

1 

0 

0 

vi 

X  = 

•V  — 

0 

0 

0 

z  ~ 

Z1  = 

Final 

X 

=  x2  = 

0 

rr 

=  z^  = 

~r 

2 

2 

0 

0 

X 

=  xn  = 

2 

2 

0 

0 

0 

z 

—  rr  — 

—  — 

2 

b  *  2  jyiA. 

These  are  shown  in  Eig.8* 
(for  small  perturbations) 


c artes ian_  perturbations 

At  the  start  of  transfer  the  relationships  are 


Sr1  =  6z^  (B.1) 

Sfc,  =  *”  6x1  (B.2) 
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8V1  =  oS1  (B.3) 

66  =  J-  6§,  +  J-  6Xl  .  (BA) 

The  above  relations  are  derived  from  geometrical  considerations  together  with 
the  assumption  that  sin  80^  s  88_j  and  oos  60^  c*  1 • 

At  the  end  of  transfer  the  inverse  expressions  can  be  derived  by  similar 
considerations : - 


6*2  - 

-r2 

(B.5) 

8s2  = 

-8r2 

(B.6) 

6%  - 

-6Y2 

(B.7) 

5l2  - 

"V2  662  +  V2  **2  • 

(B.8) 

B.3  Basic  expressions 


esian  and  -polar  derivatives 

In  polar  co-ordinates,  initial  and  final  perturbations  are  related  by  the 
following  expressions 


dr. 


dr 


dr 


dr. 


8r„  = 


6r 


dr,,  1  80. 


.  +  ~  8V  +  — i.  80 

1  80,  r1  dV,  1  80,  1 


(B.9) 


d<f>n  80  d  0 

*  R"  '  ^  CJ  '  *  ■  — i-  60 

ae1  i 


6(S2  =  ^5ri  +  wf  6vi  * 


(B.10) 


8V„ 


dY  dV  dV  dV 

6r,  +  3- ~  60,  +  6V,  +  ™  60, 

d^  1  d01  y1  dV1  1  d0_j  1 


(B.11) 


80„  = 


80  80  86  80 

5^  S  +  601  +  W1  +  667  661 


(B.12) 


Substituting  in  equations  (B.9)  to  (B.12)  for  the  expressions  of  6r^ ,  89^, 
8V^ ,  80 and  Sr^,  S02,  88^  obtained  in  equations  (B.1 )  to  (B.8)  we  can 

now  relate  start  and  end  perturbations  in  the  cartesian  axis  system,  i.e. 
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30?  o  d^2  3^2  /I  o  1 

8x_  =  -  r  v=t”  8x,  -  r_  -r— -  8z  -  r  -r^r*  I  TT"  Sz.  +  —  6x. 
2  2  dVj  1  2  1  2  dG1  \V^  1  r1  1 


r2  d#2 


r,  30,  1 


”  Sx 
1  1 

....  (B.13) 


!Tl  6o  _  r2  6  ^2  /I  60  +  1,  5  \  _  j~  2  c 

3V_j  1  3r1  1  361  ^  1  V  r1  d^1  1 


(B.14) 


av  av„  av„  .  .  \  av.  . 

2  ro  2c  2/1  cO  1  c  \  21c 

ay,  i  ar  i  ao  \v  i  r  iy  a0  r1  i 


(B.15) 


30o  30 


30o  ao 


2  _ _ j| 

.  Or. 


'2  |_ae1  ae1 


4;  •«,  *  *  *,  J  *  *,  ft 


£  U-6X, 

a1J  r1  1 


.  (B.16) 


But  by  definition 


3x2 

6x 

,  ^2 

5z. 

♦  2a 

8x. 

+  !£ 

dxl 

1 

‘  dz1 

1 

aSi 

1 

dzi 

8x 

oa 

+  lfs 

8x 

+  .ti 

ax, 

i 

1 

dzi 

1 

sSi 

1 

■*0 

8z1 

By  equating  terms  we  can  first  relate  cartesian  and  polar  derivatives  in  the 
following  dimensional  form:- 


3x? 

I 

L** 

m2  >t2- 

ri 

vei” + 

ax2 

d<p2 

8zi 

=  ~  r2 

FrT 
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iz2  =  v  r~  - — 

dz.  2  dr.  dr. 

1  Li  1 


r%  dQ 
1  2  -  avj- 


v  ra8„  ae  . 

2(2  2 


vi  td6i  a6i 


B .4  Normalised  cartesian  derivatives 

These  cartesian  derivatives  can  be  normalised,  using  the  previously 
established  normalised  polar  derivatives  and  the  various  relationships  between 
V1 ,  Vg,  r^  and  r ^  (see  Appendix  A).  The  results  are  as  follows:** 

dx2 

«*-=■  =  n+2 


221  il+sl! 

4  Vn 


ixj. 

r . 

1 


J2(l+n  )' 


Jl  1*2 

2?  . 

i  az„ 


/  2(1 +n)' 


-n  (n+2) 


3 
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V 

°1 

r1 

ax1 

V 

°1 

i!2 

r1 

ai1 

r1 

eS2 

Y 

°1 

dx1 

ri 

a§2 

V 

°1 

dz1 

^2 

ax. 

1 

^2 

dz1 

r1 

dz2 

V 

°1 

ax. 

1 

dZ2 

°1 

aZ1 

|  2n(l+n) 


3 


0 


0 


0 


2(l+n) 


n 


=  .  -25-  (1*n)2 

AnVn 
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DERIVATIVES  RELATING-  ERRORS  IN  FINAL  (NOMINALLY  CIRCULAR)  ORBIT 
TO  ERRORS  AT  START  OF  HOHMAMN  TRANSFER  -  FINAL  IMPULSE 
ASSUMED  TO  BE  APPLIED  WITHOUT  ERROR 


C  •  1  Errorsjn  the  orbital  parameters  of  a  nomir 

in  terms  of  injection  errors 


We  shall  consider  only  the  effects  upon  two  orbital  parameters,  namely 
mean  radius  and  eccentricity* 


General  expressions  for  these  parameters,  in  terms  of  cut-off  conditions 
have  been  given  in  many  forms,  but  equations  in  a  convenient  manner  for  our 
present  purposes  can  be  found  in  Refs. 5,  6*  They  are: 

Mean  radius 


Eccentricity 


where 


(l-p)  +  p(2~p)  sin  6 


(C.1) 


(C.2) 


(C.3) 


is  a  non-dimensional  parameter  already  introduced  in  Appendix  A.  r,  V  and  0 
describe  position  and  velocity  at  injection,  in  a  polar  axis  system  (see  Fig.l). 

For  injection  into  a  perfectly  circular  orbit  p=  1  and  0=0.  Let  us  now 
consider  the  effects  of  small  errors  in  r,  V  and  0  upon  the  resulting  errors  in 
a  and  e* 

Let  the  small  perturbations  be  6r,  5V,  60  respectively.  Then,  from 
equation  (C#3)j  and  neglecting  second  order  terms, 


(c. 

where  Vq  =  Vp/r  is  circular  orbital  speed  at  radius  r* 

Substituting  in  equation  (C.l)  we  get 


„  6r  26V 
1  +  —  + 
r  V 

o 
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a  +  6a  =  r  +  6a  = 


1  - 


r  + 
6r 
r 


6r 


26V 
V 


Hence 


6a 


=  26r  + 


26  Vr 
V 

o 


(c.5) 


5e 


(C.6) 


It  will  He  noted,  that  6a  is  independent  of  60.  This  is  to  be  expected  since 
a  is  a  direct  measure  of  the  total  energy  of  the  orbit,  which  depends  upon  V 
and  r,  but  not  upon  6,  the  direction  of  the  velocity  vector. 

In  equation  (C.3),  the  positive  sign  is  taken  outside  the  square  root. 
This  is  because  injection  perturbations  can  only  produce  positive  eccentricity; 
negative  eccentricity  has  no  physical  significance. 


C.2  Relation  between  perturbations  a.t  end  of  Hohmann  transfer  and  those  at 
injection  into  final  orbit  ” 

We  have  already  seen,  in  Appendix  A,  that  the  conditions  at  the  end  of 
the  Hohmann  transfer  are  as  follows : - 


Quantity 

Nominal  value 

Perturbation 

V 

V2 

W2 

0 

0 

se2 

r 

6r 

. .  2 

<t>  !  * 

! 

Following  the  application  of  an  impulsive  velocity  at  the  end  of  transfer 
the  position  parameters  (r,  <f)  will  not  change,  but  those  describing  velocity 
(V,  0)  will,  and  the  resulting  velocity  vector  must  be  defined  before  we  can 
formulate  the  injection  errors  into  the  final  (nominally  circular)  orbit.  The 
vector  addition  of  velocities  at  the  end  of  transfer  is  described  in  Fig, 9. 
Vector  1  represents  the  velocity  at  end  of  transfer  (V  =  Vg  +  6V2,  0  =  6©2). 
Vector  2  represents  the  impulse  added,  and  since  this  is  assumed  to  be  perfect 

(i.e.  the  same  as  that  required  for  nominal  transfer)  its  magnitude  is  (V  **V  )• 

°2  2 
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However,  its  direction  can  be  defined  in  two  ways,  both  of  which  are  correct  for 
the  nominal  transfer  case*  The  first  way  is  to  define  the  direction  of  impulse 
relative  to  the  local  horizontal  at  the  end  of  transfer*  Obviously  the  impulse 
(Vo  -V^)  should  be  applied  parallel  to  the  local  horizontal  and  we  shall  assume 

that  in  the  perturbed  transfer  it  is  still  applied  parallel  to  the  local  horizon¬ 
tal  (such  would  be  the  case  if,  in  practice,  the  apogee  rocket  were  aligned  by 
means  of  horizon  sensors)*  The  second  way  is  to  define  the  direction  of  the 
impulse  as  fixed  in  spaoe,  and  to  assume  that  this  inertial  direction  is  main¬ 
tained  even  in  the  perturbed  transfer  condition.  (This  would  correspond  to  the 
practical  case  in  which  an  inertial  reference,  such  as  a  gyro  or  spin  stabilisa¬ 
tion  was  used  to  align  the  apogee  rocket.)  Clearly  the  space-fixed  direction 
along  which  the  impulse  (Vq  -V2)  would  act  is  parallel  to  the  local  horizontal 

at  the  end  of  nominal  transfer,  but  in  general  this  inertial  direction  does  not 
coincide  with  the  local  horizontal  in  the  perturbed  transfer  condition. 

Referring  to  Fig# 9  we  see  that  the  velocity  vector  after  adding  the  impulse 
(vector  3,  which  is  the  sum  of  vectors  *1  and  2)  has  the  following  properties 

Magnitude 

(V  +  8V  )  +  (V  -  V  )  =  V  +  8V 

2  2  o2  2  o2  2 

(assuming  S02,  are  smaH  so  that  cos  60^,  cos  8^2  =  l). 

Vertical  component  (relative  to  local  horizontal) 


(v2  -  6V2)  502  ♦  K(V  -  V2) 


where  K  =  0  corresponds  to  the  case  in  which  the  impulse  is  aligned  along  the 
local  horizontal 

and  K  =  1  corresponds  to  the  case  in  which  the  impulse  is  aligned  in  a  space- 
fixed  direction. 

The  resulting  climb  angle  50^  is  given  by 


80^ 


vertical  velocity  component 


53 
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X*  c » 


%  6Q2  +  K  (iWp^)  6^2 


(neglecting  second  order  terms)* 

We  can  now  tabulate  all  four  perturbations  at  injection  into  the  final  orbit;- 


Quantity 

Nominal  value 

Perturbation 

V 

CvJ 

w2 

e 

0 

S©2  =  ^2  592  + 

r 

r2 

6r2 

i  * 

I 

S02 

C  »3  P>.elations  between  errors  at  start  of  transfer  and  errors  in  final  orbit 
parameters 


C.3.1  Mean  radius 

Using  equation  (C . 5 )  and  the  tabulation  of  perturbations  at  injection 
into  the  final  orbit  we  can  derive  the  mean  radius  error  of  the  final  orbit, 

i.e. 


6a  =  26r2  +  6V2  .  (C*7) 

°2 

Consider  as  an  example,  the  effect  of  a  perturbation  6V^  at  the  start  of  transfer. 
Then 


8r„ 


^2 

dV. 


6V. 


6V„ 


dV 

— £  8V 

dVj  1 


where  the  derivatives  d^/dV^ ,  dV^/dY  have  been  derived  in  Appendix  A. 
Substituting  these  expressions  in  equation  (C . 7)  we  get: 
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/  dr  2r  aVv 

6a  .  jL  6V  .  (2  J  +  -v-2-  3^)  W 

1  \  1  o2  1/ 


so  that 


>  dr  2r0  d¥ 

_da  2  2  2 

ev7  d¥.  V  d¥ 

1  1  o2  1 


(C.8) 


This  relationship  between  derivatives  can  be  expressed  in  a  more  convenient 
non-dimensional  form 


rue, 


V 

fi 

da 

-  2  ( 

¥ 

-°1 

av 

+  nVn 

&)1 

r1 

avi 

"  V 

Ai 

avd 

V>  V  J 

(C.9) 


This  expression  is  derived  by  remembering  that 

2 


V 


'1 


V 


=  n  (see  Appendix  A)  • 


Similar  expressions  relating  mean  radius  error  to  other  perturbations  at  the 
start  of  transfer  can  also  be  derived  by  a  similar  process*  The  results  are 
as  follows: 


(c.10) 


6a 

dr 


=  2 


1 


dr?\  /ri  dY 9 

—  ]  +  ni/n  ( sr—  — ~ 

V  4  0<1  1/  J 


(C.11) 


1  r1 


da 


rcr-  =  2 


I  ff2 

7  ^ 


nVn 


.  dV 

1  mjz 

V  80 

°i  1 


(c.12) 


Two  points  are  of  immediate  interest*  Firstly,  mean  radius  error  is,  to  a  first 
order,  independent  of  K  (the  method  of  aligning  the  impulse)  since  the  direction 
of  this  impulse  only  effects  vertical  velocity  at  final  injection.  Secondly 
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0  for  all  values  of  n 


1' 

are  all  zero. 

The  results  of  equations  (C .9)  and  (C.11-)  are  plotted  in  Figs. 14  and  19 
respectively. 

C.3.2  Eccentricity 

Using  equation (C. 6)  and  the  tabulation  of  perturbations  at  injection 
into  the  final  orbit  we  have  the  eccentricity  of  the  final  orbit,  i.e. 


since  in  Appendix  A  it  was  shown  that  dr^/dQ^ ,  av^/ 


and  dVg/d^ 


ir  26V  v  2 


6e  ss 


..  \r2 


O  V  \  c.  /  KC 

7  +  ~j  +  (^2  862  +  K  bWF2]  ^2) 


(C.13) 


Consider  now  the  effect  of  a  perturbation  6Vj  i- 


we  have 


6r, 


3r2 

— £  5V 

av,  1 
1 


av_ 


6V„  = 


1i?7sv1 


60  „  = 


d02 

8V 

av.  1 
1 


and 


6?4 


*0, 

— -  6V 

av1  ovi 


G-iving 


6e  =  ±6V 


1 


N 


^2 

2  V2 


!!s 

av. 


\ 2  /  30  3 \  2 

7)  -  if  *  K  t^2]  if)  •  <*•*> 


The  sign  on  the  right  hand  side  of  equation  (C.14)  needs  some  explanation. 
Clearly,  since  negative  eccentricities  cannot  exist,  the  positive  sign  must 
apply  if  6V,j  is  positive  and  the  negative  sign  is  used  if  6V^  is  negative. 
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provided  that  the  eccentricity  is  obtained  from  this  derivative  using 


rather  than 


As  in  the  case  of  mean  radius  error,  the  eccentricity  derivative  can  be 
expressed  in  the  dimensionless  form 


...  (C.15) 


and  similarly 


...  (C.1 6) 


...  (C.18) 
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Equations  (C.15)  to  (C.l8)  are  plotted  in  Figs. 15,  17 t  20  and  22.  It  should  be 
noted  that  (de/d^ )  for  the  case  of  the  apogee  impulse  aligned  to  the  local 
horizontal  (K=  0)  is  zero,  since  the  only  non-zero  term  in  equation  (C.l8)  is 
(see  Appendix  A), 


C.4 


Impulsive  velocities  required  to  null  mean  radius 
errors  of  final  orbit 


C.4.1  Nulling  radius  error 

Mean  radius  error  may  be  nulled  by  application  of  a  post  injection 
tangential  impulse.  For  an  orbit  of  small  eccentricity  it  is  of  no  consequence 
at  which  point  on  the  orbit  the  impulse  is  applied,  and  provided  it  is  applied 
along  the  direction  of  motion  the  resulting  change  in  mean  radius  is  given  by 

8a*  =  .  (C.19) 

V  a 
°2 

This  formula  is  provided  by  an  examination  of  equation  (C.5)#  6u  is  taken 

a 

here  as  the  corrective  tangential  velocity  and  in  order  to  null  any  initial 
error  6a,  5a1  must  be  made  equal  and  opposite  to  6a* 


Therefore 


6u 

a 


(C.20) 


6a  appears  in  the  modulus  sign  since,  for  practical  purposes,  we  are  at  present 

concerned  with  evaluating  the  magnitude  of  the  impulsive  correction  rather  than 

its  vectorial  value*  Of  course,  if  the  initial  error  6e  is  positive  6u  must 

a 

be  applied  in  a  retrograde  manner,  while  if  it  is  negative  the  required  correc¬ 
tion  will  be  forward. 

Using  6a  =  (da/dV^)  6V^ ,  for  an  initial  velocity  perturbation  and  substitut¬ 
ing  this  in  equation  (C.20)  we  get 


so  that 


where 


6u 

a 


It  K 
1 


8u 

a 


du 

a 


8u 

a 


(C.21) 


(C.22) 
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We  note  that  Su  must  be  evaluated  using  the  modulus  of  the  input  perturba- 
a 

tion.  Again,  this  is  because  we  are  concerned  with  evaluating  the  magnitude  of 
the  velocity  correction,  which  is  positive  for  either  positive  or  negative  inputs. 

Equation  (C.22)  may  be  expressed  in  the  following  non-dimensional  form:- 


and  similarly, 


du 


a 


(C.23) 


(C.24) 


1 

2nVn 


(C.25) 


.  0u 

(X.  J>§l\ 

l  a  1  i 

V  d6  2nVn  * 

°i  1 

Vri  dV 

(C.2 6) 


These  results  are  plotted  in  Figs. 16  and  21.  du  /<30.  and  Su  are  both  zero 

Si  l  Si  I 

for  all  values  of  n. 

C#4# 2  Mulling  eccentricity  error 

A  tangential  impulse  can  also  be  used  to  null  a  small  eccentricity  error. 
The  magnitude  of  such  an  impulse,  6u  is  given  by 

V 

°9 

6u  =  ~  6e  .  (C.27) 

e  c. 

This  result  is  seen  from  equation  (C.6).  It  is  interesting  to  note  that  if 

eccentricity  error  was  corrected  by  a  vertical  impulse,  rather  than  a  tangential 

one,  the  velocity  required  would  be  twice  as  large,  since  50  =  vertical 

velocity/V  •  Eccentricity  error  can  be  reduced  by  either  a  forward  impulse  at 
°2 

apogee  or  a  retrograde  impulse  at  perigee. 

From  equation  (C*27)  we  can  see  that 
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3u  o„ 

_ _ e  Je_ 

av,  av  *  2 

i  i 


(C  *  28) 


Again  we  must  remember  to  evaluate  Su^  from  the  relation 


rather  than 


—2.  5Y 
dV1  V1 


Equation  (C.28)  can  be  expressed  in  the  non-dimensional  form:- 


and  similar ly. 


-JL  =  4 Vv 

dV1  2vn  ^  0l  dV1 


(C.29) 


IV  30, 
I  v  o . 

L _ 1 _ 


•r,  3u  ■ 
1  _ e 

V  dr. 

‘  °1  1 


1  /de 

iTn  ( ae . 


X/r 

2Vn  I  1  dr 


(C.30) 


(C.31) 


'  j,  3u  • 
i  _ e 

V  30, 

,  o,  y\. 

1 


1  / de 

2Vn  i  30 


(c.32) 


It  should  be  noted  that  there  is  no  need  to  include  a  modulus  sign  on  the  right 
hand  side  of  equations  (C.29)  to  (C.32),  since  eccentricity  derivatives  are 
already  positive.  These  derivatives  are  plotted  in  Figs. 16,  18,  21  and  23« 

C.4.3  Simultaneous  nulling  of  radius  and  eccentricity  errors 

In  general  this  type  of  correction  will  require  two  tangential  impulses 
8u^  and  6u^  (see  Fig.10).  The  impulses  must  be  allocated  in  such  a  manner  that: 

(i)  the  algebraic  sum  (taking  into  account  whether  the  impulses  are 
forward  or  retrograde)  must  null  the  radius  error, 

(ii)  the  two  impulses  must  be  split  between  apogee  and  perigee  so  that 
the  nett  eccentricity  correction  equals  the  initial  eccentricity  error. 

TT.TdFTO  IJ23AB.1 
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In  Fig.10(a)  the  perturbed  orbit  lies  wholly  outside  the  desired  circular 
orbit  and  hence  both  impulses  are  retrograde  (or  forward).  Thus  the  total 
impulse  requirement  in  this  case  is  dictated  by  radius  correction  and  the 
eccentricity  correction  can  be  accommodated  within  this  requirement. 

In  the  other  case  (Fig*10(b))  the  perturbed  orbit  intercepts  the  desired 
one,  and  both  impulses  are  instrumental  in  nulling  eccentricity.  Thus  the 
total  requirement  is  dictated  by  the  initial  eccentricity  and  mean  radius 
correction  can  be  accommodated  wdthin  this  requirement# 

The  condition  that  the  perturbed  ellipse  does  not  intercept  the  desired  one 
is  that 


or 


i.e. 


(l~6e)  +  6a  > 


(l+6e)  +  6a  < 


6a  |  >  ^ 


(C.33) 


Comparison  of  this  equation  with  equations  (C.20)  and  (Co27)  shows  that  this  is 
equivalent  to  the  condition 


6u  >  6u 
a  e 


as  would  be  expected. 

To  summarise  we  can  suppose  that  in  general  ou  can  be  either  greater  or 

a 

less  than  Su^,  depending  upon  the  nature  of  the  perturbed  orbit.  Both  mean 

radius  and  eccentricity  errors  can  be  nulled  by  application  of  two  tangential 

impulses  whose  arithmetical  sum  equals  the  larger  of  6u  and  6u  ,  provided  the 

e  a 

impulses  are  applied  at  the  correct  p>laces. 

It  will  be  seen  by  examining  the  curves  of  Figs. 1 6,  18,  21  and  23  that 
8ue  may  indeed  be  larger  or  smaller  than  8u^,  depending  upon  the  nature  of  the 
perturbation. 
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SYMBOLS 


a 


e 

K 


T 

o 


V 

v0  =Vi7? 


X 


1 


mean  radius  of  orbit 
eccentricity 

coefficient  defining  method  of  alignment  of  apogee  impulse 
transfer  ratio 


dimensionless  parameter  relating  speed  to  circular  orbital  speed 

geocentric  radius 

orbital  period  of  transfer  orbit 

magnitude  of  post-orbital  velocity  correction  required  to  null  mean 
radius  error 

magnitude  of  post-orbital  velocity  correction  required  to  null 

eccentricity  error 

velocity  magnitude 

ciroular  orbital  speed  at  radius  r 

cartesian  position  co-ordinates 


6 

Subscripts 

1 

2 


flight  path  angle  relative  to  local  horizontal 
angular  range 

earth’s  gravitational  constant 

point  at  which  Hohmann  transfer  starts 
point  at  which  Hohmann  transfer  ends 
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FIG.  16  VELOCITIES  REQUIRED  TO  CORRECT  MEAN  RADIUS 
&  ECCENTRICITY  ERRORS  OF  FINAL  ORBIT,  DUE  TO 

VELOCITY  ERROR  Sv,  AT  START  OF  TRANSFER 


r-ooose 


SPA/P/IG45 


Fig. 17 


£3456789  lO  VALUE  FOR 


h2  (nautical miles!  for  h,-30onm 


FIG.I7  ECCENTRICITY  OF  FINAL  ORBIT  DUE  TO  CLIMB 
ANGLE  ERROR  6  ©i  AT  START  OF  TRANSFER 
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